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a s Determinations from Jets and Scaling Violations at HERA 

T. Kluge 

DESY, Notkestr. 85, 22607 Hamburg, Germany 

A review is given on recent a s determinations from the HI and ZEUS Collaborations. These are 
based on measurements of jet cross sections, event shape variables, as well as on the observed scaling 
violation of the structure function F%. A HERA average on a s {mz) is presented, in comparison 
with world mean values. 
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The strong coupling constant a s is the single free 
parameter of QCD, the knowledge of its value is es- 
sential when predicting virtually every cross section 
for high energy collisions of elementary particles. E.g. 
multi jet states pose the standard model background 
for various searches for new physics at the LHC. Thus 
it is essential to know a s to as high precision as possi- 
ble. The determination of the strong coupling is more 
difficult compared to other elementary forces due to 
the confinement, i.e. one cannot observe the carriers 
of color charge, the partons, directly like electrons in 
QED. Before the startup of LEP and HERA ajjreci- 
sion of only around the 10% level was reached 

a s (m z ) = 0.11 ±0.01. 
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Much progress has been made since then, thanks to 
experimental data of higher precision, larger range in 
scale and the inclusion of a variety of processes, such 
as e + e~ annihilation, hadron-hadron collisions, ep 
scattering, heavy particle decays as well as advances 
in the theory. 

HERA, the only accelerator for e ± p collisions at 
high center-of-mass energies y/s — 320 GeV, is recog- 
nised as today's precision tool for QCD investigations. 
With the HI and ZEUS experiments two general pur- 
pose detectors explore the manifold aspects of QCD. 
The role of HERA for determinations of the strong 
coupling constant lies in precision as well as in com- 
plementarity to other environments. In general two 
effects are exploited for this: scaling violations of 
structure functions (inclusive) and hadronic final state 
characteristics (exclusive) . 

The cross section of the neutral current (NC) in- 
teraction ep — > eX (Fig. ^) is defined in terms of 
three kinematic variables Q 2 (the photon virtuality), 
Bjorken x and y, where y quantifies the inelasticity 
of the interaction. The kinematic variables are re- 
lated via Q 2 = sxy, where s is the ep centre-of-mass 
energy squared. Over most of the large kinematic 
domain at HERA the dominant contribution to the 
cross section comes from pure photon exchange as 
expressed by the structure function F%(x, Q ), which 
is in the quark parton model directly related to the 
sum of the quark and antiquark densities in the pro- 
ton, i.e. xF 2 = J2i e l(<li( x >Q 2 ) +Qi(x,Q 2 ))- In the 
quark parton model F2 is independent of Q 2 , a fea- 
ture known as scaling, consequently the observation 
of scaling violations is evidence for QCD. The accu- 




FIG. 1: The Born contribution to neutral current ep scat- 
tering. 



racy and kinematic coverage of the HERA neutral and 
charged current cross section data 0- 03- [3] - provide for 
a clear demonstration of scaling violations, e.g. shown 
in Fig. |2 The magnitude of the scaling violations de- 
pends on a s as well as on the gluon density in the pro- 
ton. The HI Collaboration [j| has determined a s {mz) 
simultaneously with the gluon density in a QCD fit to 
the derivative of F 2 with respect to Q 2 (Fig. |3J , using 
HI and BCDMS proton data. A value of the coupling- 
constant 

a s (m z ) = 0.1150 ± 0.0017(exp.) ±£$$5 (mod.) 

± 0.0050(th.) 

is obtained. The theory error, estimated by varying 
the renormalisation scale and the factorisation scale 
by a factor of four, is found to be significantly larger 
than the experimental uncertainty. This uncertainty 
is expected to be reduced significantly in NNLO per- 
turbation theory. The value obtained for a s (mz) is 
nearly independent of the chosen parameterisation for 
the parton densities. 

A more direct sensitivity to a s is achieved when in- 
cluding the hadronic final state of ep scattering. An 
example is the cross section for events where the final 
state contains more than one hadronic jet (not count- 
ing the remnant from proton dissociation) , which van- 
ishes in absence of QCD effects (Fig. 0J). In dccp- 
inclastic scattering (DIS, Q 2 > lGeV 2 ) the jets are 
usually reconstructed in the Breit frame of reference, 
where by definition no recoil from the scattered beam 
electron occurs. Jet cross sections are complementary 
compared to inclusive NC cross sections in that sense 
that the analysis offers direct sensitivity to a s , on the 
other hand the measurement as well as the theoreti- 
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FIG. 2: The proton structure function F2 measured by 
the HI Collaboration compared with the corresponding 
Standard Model expectation determined from the HI PDF 
2000 fit (error bands). Also shown are the Fi data from 
BCDMS and NMC, which are not used in the fit. 
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FIG. 3: The derivative (dF2/d In Q 2 ) x plotted as functions 
of x for fixed Q 2 , for the HI data (symbols) and the QCD 
fit to the HI data, for Q 2 > 3.5 GeV 2 (solid lines). The 
error bands represent the model uncertainty of the QCD 
analysis. 



FIG. 4: The lowest order diagrams for two jet production 
in the Breit frame of reference. 
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FIG. 5: Measured inclusive jet cross section da/dE t com- 
pared with NLOJETH — h predictions corrected for hadro- 
nisation effects. The bands show the theory uncertainty 
by varying the renormalisation and factorisation scale by 
a factor of two. The statistical uncertainty is shown as 
the inner error bar, while the total error is denoted by the 
outer error bar. 



cal treatment is more involved, resulting to in general 
larger uncertainties of the data points. The HI and 
ZEUS Collaborations published results of inclusive jet 
and dijet cross sections [(JlI- In a recent analysis the 
HI Collaboration presented a measurement of the in- 
clusive jet cross section in DIS j| where the jets are 
defined by the k t algorithm in the Breit frame of refer- 
ence. The phase space spans high values of the photon 
virtuality 150 < Q 2 < 5000 GeV 2 . In Fig. the single 
differential cross section as function of the transverse 
energy E t of the jets is shown. A NLO calculation, 
corrected for hadronisation effects, gives a good de- 
scription of the data over the full E t and Q 2 range. 
Hadronisation corrections are determined with Monte 
Carlo event generators, they are found to be typically 
around 10% from unity. A fit of the strong coupling 
constant was performed by calculating the theory for 
several values of a s (mz) and interpolating the cross 
section in between. This interpolation allows for a 
mapping of the measured cross sections together with 
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FIG. 6: Results for a s (mz) and a s (Et), where the a s (E t ) 
values have been evolved from the measured a s (mz) val- 
ues for each bin in E± using the two loop solution of the 
renormalisation group equation. The solid curve shows the 
result of evolving the averaged a 3 (mz), while the dashed 
curves indicate the respective error band. The horizontal 
line with the shaded band corresponds to a world average. 



their uncertainties onto a corresponding a s (mz) in- 
terval. The result are shown in Fig. as a function 
of the scale E t . The running of the strong coupling is 
clearly visible, and the values are found to be consis- 
tent within errors to a world mean value. Averaging 
the four individual values and taking into account cor- 
related systematical errors yields 

a s (m z ) = 0.1197 ± 0.0016(exp.) 



+0.0046 
-0.0048 



(th.). 



Again the error from theory is significantly larger than 
the experimental uncertainties. The largest contribu- 
tion to the latter comes from the hadronic energy scale 
of the main calorimeter. 

For inclusive jet cross sections every single jet is 
counted, alternatively the frequency of multi jet events 
allows for another handle on the size of a s . One possi- 
bility is to study the cross section ratio of three-jet and 
two-jet events, which is approximately proportional to 
a s . By construction a partial cancellation of errors, 
e.g. from the uncertainty of the gluon density in the 
proton, is expected. There are results from HI as 
well as from ZEUS available, which investigate 
the multi jet ratio ratio in DIS. The latter is shown as 
a function of the photon virtuality Q 2 in Fig. [7| The 
data are compared to a calculation at NLO precision, 
which has been corrected for hadronisation effects. A 
good description within errors is obtained of both, the 
shape and the normalisation of the measured ratio. 
The total experimental and theoretical uncertainties 
are about 5% and 7%, respectively. These uncertain- 
ties are substantially reduced with respect to those of 
the individual di- and trijet cross sections. In par- 
ticular, at lower Q 2 the theoretical uncertainties are 
reduced by as much as a factor of four. In the upper 
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FIG. 7: (a) The ratio of inclusive trijet to dijet cross sec- 
tions as a function of Q 2 . The predictions of perturbative 
QCD in next-to-leading order are compared to the data, 
(b) shows the ratio of the data to the predictions. The in- 
ner error bars represent the statistical uncertainties. The 
outer error bars represent the quadratic sum of statisti- 
cal and systematic uncertainties not associated with the 
calorimeter energy scale. The shaded band indicates the 
calorimeter energy scale uncertainty. 
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FIG. 8: (a) The ratio of inclusive trijet to dijet cross sec- 
tions as a function of Q 2 . The predictions of perturbative 
QCD in next-to-leading order using five sets of CTEQ4 
PDF are compared to the data, (b) shows the a s (Mz) 
values determined from the ratio of inclusive trijet to dijet 
cross sections in different regions of Q 2 . The shaded band 
indicates a world average value of a 3 (Mz)- The dashed 
error bars display the theoretical uncertainties. 
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diagram of Fig.[5]the sensitivity to the strong coupling 
is highlighted. The theory prediction is shown for dif- 
ferent values of a s (mz)- A method similar to what 
was used for the inclusive jets yields values of a s (mz) 
as shown in the lower diagram. The mean value 
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a s {m z ) 



0.1179 ±0.0013(stat.) lo.oo46( ex P-) 

+0.0064/,, s 
-0.0046 



is found to be compatible with what was extracted 
from the inclusive jets, where the latter exhibit a 
somewhat smaller uncertainty. On the other hand, 
the ratio of multi jet cross sections allows for an exten- 
sion of the phasespace to lower Q 2 (10 GeV 2 compared 
to 150 GeV 2 ), due to cancellations of the theoretical 
uncertainties. 

In the regime of photo production at even lower 
Q 2 < 1 GeV 2 an additional resolved component to 
the cross section needs to be taken into account. In 
that case the photon acts like a hadron as a source of 
partons, which in the theory is reflected by the intro- 
duction of a second set of parton density functions. 
The large cross section for photo production allows to 
extend the measurement to higher jet transverse en- 
ergies. An analysis by the ZEUS Collaboration 11] 
spans 17 < E t < 95 GeV, shown in Fig. The mea- 
sured da I dE^ falls by over five orders of magnitude 
m this E^ range. Fixed-order QCD calculations are 
compared to the data, where the LO QCD calculation 
underestimates the measured cross section by about 
50% for E^ < 45 GeV. The calculation that includes 
NLO corrections gives a good description of the data 
within the experimental and theoretical uncertainties 
over the complete E 3 ^ range studied. 

A value of a s (mz) is obtained by combining all the 



E^ regions 

a s (m z ) = 0.1224 ± 0.0001(stat. 



+0-0022/ s 
-0.0019 K^P-J 

+0.0054/,, x 
-0.0042 



which is consistent with the determinations from jet 
production in NC DIS. The largest contribution to 
the experimental uncertainty comes from the jet en- 
ergy scale and amounts to ±1.5% on a s (mz)- Among 
the contributions to the theoretical uncertainty on the 
strong coupling, the largest one is due to terms be- 
yond NLO, which is estimated as ^3' 2 % by varying the 
renormalisation and factorisation scales by a factor of 
four. Uncertainties connected to the proton and pho- 
ton PDFs and to the hadronisation were significantly 
smaller. 

Event shape variables play an important role in 
QCD studies of the hadronic final state. They are 
a class of topological observables which are defined 
as positive real numbers calculated from the four vec- 
tors of all hadronic final state part icles. Both, the HI 
and ZEUS Collaborations [l2L Il3j| published recently 
new investigations of event shape variables, motivated 
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FIG. 9: The inclusive jet cross section da/dE 1 ^ in pho- 
toproduction. The uncertainty associated with the abso- 
lute energy scale of the jets is shown as a shaded band. 
The LO (dashed line) and NLO (solid line) QCD parton- 
level calculations corrected for hadronisation effects are 
also shown. The lower diagram shows the fractional differ- 
ence between the measured da/dE 1 ^ and the NLO QCD 
calculation; the hatched band shows the uncertainty of the 
calculation. 



by new developments from the theory side. There 
are several differences compared to jet cross sections: 
larger statistics due to the semi-inclusive nature of 
event shapes, reduced experimental systematic uncer- 
tainties from hadronic energy scales and larger hadro- 
nisation effects. The hadronisation needs a special 
treatment, in recent DIS studies by analytic power 
corrections. In general event shapes variables are de- 
fined to be equal to zero without QCD effects and are 
of order one for multi jet configurations, which can 
occur due to multiple gluon radiation. For differential 
distributions, fixed order calculations alone proved to 
be not sufficient, hence soft gluon resummations are 
supplementing the NLO predictions. The HI Collab- 
oration [l2T ] published distributions of five event shape 
variables, shown in Fig. ^| Towards higher scales Q 
the distributions are more peaked near zero (the value 
for the quark parton model) , which gives evidence for 
asymptotic freedom of QCD. This distributions are 
compared to calculations based on fixed order and 
matched resummed parts. To take into account hadro- 
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FIG. 11: Fit results to the differential distributions of r, B, 
po, re and the C-parameter in the (a s ,ao) plane. The la 
contours correspond to \ 2 = Xmin + 1; including statistical 
and experimental systematic uncertainties. The value of 
a s (vertical line) and its uncertainty (shaded band) are 
taken from |l4j|. 



FIG. 10: Normalised event shape distributions corrected 
to the hadron level for the to, t and B variables. The 
measurements are compared with fits based on a NLO 
QCD calculation including resummation (NLL) and sup- 
plemented by power corrections (PC). The fit results are 
shown as solid lines and are extended as dashed lines to 
those data points which are not included in the QCD fit. 



nisation Dokshitzer/Webber power corrections (PC) 
have been used, which depends on the parameter a 
representing an effective strong coupling constant in 
the infrared regime. An overall good description is 
obtained for part of the phase space (higher Q and 
moderate event shape values), where the theory is 
expected to be valid. Simultaneous fits of a s (mz) 
and the power correction parameter ao are shown in 
Fig. ^] An average value of 

a s (m z ) = 0.1198 ± 0.0013(exp.) 



+0.0056 
-0.0043 



(th.) 



is obtained, which is consistent with the results from 
jet and inclusive DIS cross sections. The fit was also 
performed separately for all scales covered by the data, 
see Fig. El where the asymptotic freedom of QCD 
is clearly demonstrated. Due to the more inclusive 
definition compared to jets, a larger range in scale is 
accessible for the event shape analysis. 

Since inclusive DIS and jet analyses offer different 
sensitivity to the PDFs of the proton and a s , it is 
desirable to have a combined QCD analysis based on 
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FIG. 12: The strong coupling a s as a function of the scale 
Q from an average of the results obtained by fitting the 
differential event shape distributions. The errors represent 
the total experimental uncertainties. A value of a s (mz) is 
indicated in the plot, determined from a fit to the a s (Q) re- 
sults using the QCD renormalisation group equation. The 
fit curve is shown as the full line. The inner (outer) shaded 
band represents the uncertainty of the fitted a s (Q) from 
experimental errors (the renormalisation scale variation). 
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FIG. 13: The \ 2 profile as a function of a s (mz) for the 
ZEUS-JETS-q s fit (black dots) and for a similar fit not 
including the jet data (clear dots). The ordinate is given 
in terms of the difference between the total x 2 an d the 
minimum x* ''■> for each fit. 
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FIG. 14: ZEUS- JETS fit compared to ZEUS DIS jet data. 
Each cross section has been multiplied by the scale factor 
in brackets to aid visibility. 



both classes of observables, simultaneously fitting the 
PDFs and the strong coupling. The addition of jet 
data allows an extraction of a s (mz) that is less cor- 
related to the shape of the gluon density in the pro- 
ton. Such an investigation was presented by the ZEUS 
Collaboration |15| . A clear improvement on the un- 
certainty of a s (mz) by adding jet data to the inclu- 
sive fit is evident from Fig. El The included jet data 
comprise DIS as well as photo production data, an ex- 
ample is shown in Fig. 1141 The result of the combined 
fit for the strong coupling yields 

a s {m z ) = 0.1183 ± 0.0028(exp.) ± 0.0008(mod.) 

± 0.005(th.), 

where in contrast to the determinations given before 
the theory uncertainty was estimated by varying the 
scale by a factor of two (instead of four). 

From what was presented in this article so far it is 
evident that there is a great diversity of a s {mz) deter- 
minations from the HI and ZEUS Collaborations. An 
HERA average from selected values was constructed 
in |l6j . The individual results which entered the av- 
erage are listed in Fig. ED All of these are consistent 
with each other and the theory uncertainty always 
dominates the total error. The average was build by 
taking into account the known correlations in each 
experiment and assuming conservatively that the the- 
oretical uncertainties arising from the terms beyond 
NLO arc fully correlated. The HERA average 

a s {m z ) = 0.1186 ± 0.0011(exp.) 

± 0.0050(th.) 



th. uncert. . 
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FIG. 15: Summary of a s (mz) determinations at HERA 
compared with a world average. 
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FIG. 16: Summary of measurements of a a {Q) as a function 
of the respective energy scale Q |l7j. Open symbols indi- 
cate (resummed) NLO, and filled symbols NNLO QCD 
calculations used in the respective analysis. The curves 
are the QCD predictions for the combined world average 
value of a 3 (mz), in 4-loop approximation and using 3- 
loop threshold matching at the heavy quark pole masses 
M c = 1.5 GeV and M b = 4.7 GeV. 



is compatible with world means (e.g. from and 
offers competitive precision. 

This HERA average was included in a more recent 
world average |17J > 

a s {m z ) =0.1189 ±0.0010, 

where for the first time jet data from HERA were con- 
sidered. The individual contributions as a function of 
their scale are shown in Fig. ^| The total error on 
this average is quoted to be below 1%, a large im- 
provement compared to previous world means, which 
could be reached by including contributions from tau 
decays and lattice QCD. 

The uncertainties of a s (mz) determinations from 
HERA data alone are somewhat larger, still compared 
to e + e~ annihilation and hadron collisions they are 
competitive. Moreover, the observed compatibility of 
the results from HERA and other processes is essential 
to support the universality of QCD. 
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